Stable Fulde-Ferrell-Larkin-Ovchinnikov pairing states in 2D and 3D optical lattices 
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We present the study of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) pairing states in the p- 
orbital bands in both two and three-dimensional optical lattices. Due to the quasi one-dimensional 
band structure which arises from the unidirectional hopping of the orthogonal p-orbitals, the pairing 
phase space is not affected by spin imbalance. Furthermore, interactions build up high dimensional 
phase coherence which stabilizes the FFLO states in 2D and 3D optical lattices in a large parameter 
regime in the phase diagram. These FFLO phases are stable with imposing the inhomogeneous 
trapping potential. Their entropies are comparable to the normal states at finite temperatures. 
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The FFLO phases are a class of exotic Cooper pairing 
states exhibiting non-zero center of mass momenta ■> 
which occur in spin imbalanced systems with mismatched 
Fermi surfaces. However, such states are difficult to re- 
alize in solid state systems. The strong orbital effects 
of external magnetic fields often suppress Cooper pairing 
before sizable spin polarizations are reached. Moreover, 
because only small fractions of the mismatched Fermi 
surfaces can participate pairing, the FFLO states are usu- 
ally fragile in 2D and 3D systems. In spite of indirect ev- 
idence in various heavy fermion compounds and organic 
superconductors (e.g. CeCoIn 5 [5] and A-(BETS) 2 FeCl 4 
Q), the FFLO states remain elusive. 

In the cold atom community, the search for the FFLO 
pairing states has been attracting considerable interest 
[7H27| . Spin imbalanced two-component fermion systems 
have been prepared free of the orbital effects of magnetic 
fields. However, the problem of the limited pairing phase 
space remains, thus phase separations are observed ex- 
perimentally instead of the FFLO pairing in 3D traps 
[2 a |28|. This difficulty is avoided in ID systems whose 
Fermi surfaces are points, thus spin imbalance does not 
affect the pairing phase space. Considerable progress has 
been made in quasi- ID systems of coupled optical tubes 
in Hulet's group Q, in which the partially polarized cen- 
tral regions in the tubes are observed in agreement with 
the prediction of the Bethe ansatz solution. However, 
due to the intrinsic strong quantum fluctuations in ID, 
the pairing density waves, which are the smoking gun evi- 
dence for the FFLO states, cannot be long-range-ordered 
and thus difficult to observe. 

On the other hand, orbital physics with cold atoms in 
optical lattices has received considerable attention, which 
gives rise to a variety of new states of matter with both 
cold bosons and fermions [291-134} . In particular, it has 
been recently shown that the p x ^-orbital band in the 
honeycomb lattice exhibits different properties from its 
p^-orbital counterpart of graphene. These include the 
strong correlation effects in the flat bands (e.g. Wigner 
crystallization [35| and ferromagnetism [36|), quantum 



anomalous Hall states |37l| . and the heavily frustrated 
orbital exchange physics [38|, [39| . 

In this article, we combine the realization of the FFLO 
states and the study of orbital physics with cold atoms 
together. The FFLO states can be stabilized in the 
p-orbital bands in both 2D square and 3D cubic opti- 
cal lattices. Different from the met ast able p-orbital bo- 
son systems {2jJ the p-orbital systems filled with 
fermions with the fully filled s-band are stable due to 
Pauli's exclusion principle. This work is a natural high 
dimensional generalization of the current experiments in 
Hulet's group Q- The p x (p y ,p z )-orh\tdl bands behave 
like orthogonally-crossed quasi- ID arrays due to their 
highly unidirectional hoppings. The onsite negative Hub- 
bard interactions further build up high dimensional phase 
coherence over the entire lattice. It combines the advan- 
tages of the large pairing phase space of quasi- ID systems 
and the high dimensional phase coherence. 

The anisotropic p-orbital bands possess the quasi- ID 
like structures with perfect nesting at general fillings and 
spin imbalance. For simplicity, we start with the 2D case. 
The similar physics applies to the 3D cubic lattice as well. 
We present the p-band Hamiltonian as 
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where a refers to spin index; h controls spin imbalance; 

n a (f) = pt^ffiPxAr) +Pj,a(r)Py,a(r) is the Particle 
number of spin a. Only the longitudinal a-bonding (ty) 
term is kept which describes the hopping between p- 
orbitals along the bond direction as depicted in Fig. [1] 
(a). t\\ is positive because of the odd parity of the p- 
orbitals. The transverse 7r-bonding term with the hop- 
ping integral t± is neglected, which describes the hopping 
between p-orbitals perpendicular to the bond direction as 
depicted in Fig. D](b). 

In spite of the 2D lattice structure, the p-orbital band 
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FIG. 1: (Color online) (a) and (b) describe the longitudi- 
nal hopping 1 1| -term and the transverse £j_-term of the the 
p-orbitals, respectively, (c) The pairing hopping term in Eq. 
([2]) locks the phases of two onsite intra-orbital pairings in the 
p x and p y -orbitals. 



structure of Eq. [T] remains quasi-ID-like as depicted in 
Fig. [2] (a). The p x {Py) -orbital band disperses along the 
#(2/)-direction, respectively, but does not along the ^(re- 
direction. The Fermi surfaces are vertical (p x ) and hori- 
zontal (jp y ) lines across the entire Brillouin zone. For the 
arbitrary filling and spin imbalance, the Fermi surfaces 
of spin up and down fermions have the perfect nesting. 
Consequentially, spin imbalance does not suppress the 
pairing phase volume. The high dimensional p-orbital 
systems have the same advantage as that in ID systems. 

The important feature of the 2D p-orbital systems for 
the FFLO states is that the onsite negative Hubbard in- 
teractions build up the 2D phase coherence. The inter- 
actions are represented in the standard two-orbital Hub- 
bard model as 



Hint — 



J2 U [ n ^(r)n xi (r) + n yt (r)n yi (r) 

r 

^j[s x (r) ■ S y (r) - ^n x (r)n y (r) 

r 

Y, A \pU(^pU(r)Pyi(r)Pyi:(r) + h.c] , (2) 



where U — g J dr\ip Px y (r)| 4 < and g is the contact in- 
teraction in the s-wave scattering approximation. J and 
A satisfy J = ^ < and A = f < The nega- 

tive [/-term gives rise the dominant intra-orbital singlet 
pairings in the p x and p^-orbitals, defined as 



A x (f) = (G\p xt (?)p xl (r)\G), 
A.(r-) = (G\p vt (f) Pvi (f)\G), 



(3) 



where \G) is the mean field pairing ground states. The J- 
term induces the inter-orbital singlet pairing between p x 
and p^-orbitals. However, because the Fermi surfaces of 
p x and p^-orbitals are orthogonal, the inter-orbital pair- 
ing is unfavorable. 

The pair hopping A-term in Eq. [2]can be considered as 
the internal Josephson coupling to lock the phases of two 
intra-orbital pairings A^ and A^. As a result, the motion 
of Cooper pairs are 2D-like in spite of the quasi ID-like 
single fermion hopping. To clarify the pairing symmetry, 
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FIG. 2: (Color online) (a)The nesting of the p-orbital Fermi 
surfaces ensures that all of the Fermi surfaces are paired at 
a general filling and spin imbalance. The Fermi surfaces of 
the p x (p y ) -orbit als are vertical (horizontal) lines; those of 
the majority (minority) spins are marked red (blue). Fermi 
surfaces marked with solid (dashed) lines are paired with the 
center of mass momentum ±Q, respectively. The red and blue 
arrows represents the Fermi wavevectors of spin up and down 
fermions participating Cooper pairing, (b) The 2D phase di- 
agram from the B-de G solution as chemical potential fi and 
the magnetic field h with U/t\\ = —1.5. 



we first consider two fermions on the same site to gain 
some intuition. The s-wave Feshbach resonances forbid 
spin triplet channel and induce a spin singlet pairing. In 
the spin singlet channel, their orbital wavefunctions are 
symmetric as p\ -\~Py, p x ~Py and p x p y respectively. The 
first one has energy L^+A = 4C//3, while the later two are 
degenerate with energy U — A = J = 2U/3. From this 
simple analysis, we can see that the system favors pairing 
with p\ + p^ orbital symmetry (as shown in Fig. Hfc)), 
while pairing with other two symmetries are suppressed, 
which can be verified by the numerical results below. 

We have performed calculations based on the self- 
consistent Bogoliubov-de Gennes (B-deG) solution to 
study the competition among the FFLO state, the BCS 
state and the normal state as presented in Fig. [2] (b). 
To synchronize the phases of A x (f) and A 2/ (r) on each 
site, their center of mass wavevectors in the FFLO states 
have to be the same. This can be achieved by choosing 
the pair density wavevectors along the diagonal direction 
±Q defined as Q = (5kf,5kf) where 5k f = kf 1 —kf 21 and 
kf 12 are Fermi wavevectors of the majority and minority 
spins as indicated in Fig. [2] (a). By the symmetry of 
the square lattice, Q' = ±(5kf,—5kf) are another pos- 
sible choice of pair density wave vector. We consider the 
simplest Larkin-Ovchinnikov (LO) states with one pair 
of Cooper pair momenta ±Q, with the sinusoidal order 
parameter configuration as 



A,(f)=A,(f) = |A|cos(Q-f). 



(4) 



The LO state breaks both translational and the 4-fold 
lattice rotational symmetries. We have performed un- 
biased real space B-de G calculations without specifying 
the FFLO momentum in the initial conditions but rather 
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starting from a configuration with uniform pairing. The 
FFLO momentum Q in the above analysis is obtained 
when the numerical convergence is arrived. Compared 
with the phase diagram of spin-imbalanced fermions in 
the s-orbital band, the FFLO phase in our p-orbital band 
system exists in a much larger regime in the phase dia- 
gram sandwiched between the fully paired BCS phase 
and the fully polarized normal phase. 




FIG. 3: (Color online) The B-de G solution of the p-orbital 
FFLO state in a 30 x 30 lattice with a weak confining trap, (a) 
The order parameter distribution of (f) oscillates along the 
[1 Indirection, (b) The spin density distribution S z (r) peaks 
around the gap nodes. 

Next we study a more realistic situation, the effects of 
the soft confining potential to the p-orbital FFLO states, 
by performing self-consistent real space B-de G calcula- 
tions. We consider a 30 x 30 lattice with the harmonic 
trapping potential V ex = A(r/a) 2 where A/tn = 5x 10 -3 ; 
r is the distance from the trap center; a is the lattice 
constant. The real space distribution of order parameter 
A x (f) is shown in Fig|3] (a) with the parameters cho- 
sen as h/tn = 3, U/tn = —3 and \i = 0. Clearly A x 
oscillates along the [11] -direction in agreement with the 
previous analysis. To verify Eq. [4j we further calculate 
the difference between the pairing orders in different or- 
bitals. In the bulk, the relation that A x (r) = A y (f) is 
well-satisfied. The difference between A x {f) and A y (f) 
is only important at the boundary which breaks the sym- 
metry between the p x and p^-orbitals. The spin density 
distribution s z (f) = n^(f) — n^{f) is depicted in Fig. [3] 
(b). It peaks around the gap nodes, which is consistent 
with the fact that spin polarization suppresses Cooper 
pairing. 

Next we discuss the effect of the small 7r-bonding t±, 
which has been neglected above but always exists in re- 
alistic systems. The t±-term restores the 2D nature of 
the Fermi surfaces and suppresses the perfect nesting, 
therefore it is harmful to the FFLO states. Our numer- 
ical result indicates that the FFLO state remains stable 
at small values of t±. For example, with U/tn = —3, 
h/tn =3 and \i = 0, the FFLO state survives until t±/tn 
reaches 0.12. Beyond this value, it changes to the nor- 
mal state through a first order phase transition. As cal- 
culated in Ref. [30], with the optical potential depth 
Vo/E R « 15, ty is at the order of 0.1E R and t±/t\\ « 5%. 



Increasing optical potential depth further suppresses t± , 
thus there is a large parameter regime to stabilize the 
FFLO states. 
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FIG. 4: (Color online) (a) The finite temperature phase di- 
gram for the 3D p-orbital bands with U/t\\ = —2.4, and \i — 0. 
T c /t\\ — 0.84 is the critical temperature for the BCS state, (b) 
The entropy density S/ks v.s. h/tu at a finite temperature 
of T/tu = 0.4 with parameters \i — and U/t\\ — —2.4. 



The physics of the FFLO states in the p-orbital bands 
in the 3D cubic optical lattices is similar. The 3D p- 
orbital Hamiltonian is similar to Eq. [1] and Eq. [2j and 
further augmented by a new orbital p z . The pair density 
wavevectors are along the body diagonal directions, i.e., 
the ±[111] or other equivalent directions. Similar to the 
Eq.(|4]), the FFLO state in 3D cubic optical lattice is char- 
acterized by the sinusoidal order parameter configuration 
as: 

A x (r) = A y (r) = A z (f) = \A\ cos(Q • r). (5) 

where Q = (±5kf,±5kf,±5kf) are along the body- 
diagonal directions. 

In the 3D p-orbital bands, the long range ordered BCS 
and FFLO states survive at finite temperatures, and 
mean-field theory works qualitatively well. We present 
the finite temperature phase diagram of the competing 
orders at U/t\\ = -2.4 and \i = in Fig. g] (a). The 
FFLO state can also survive to finite critical tempera- 
tures at the same order of T c . We further present the 
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entropy S v.s. h for different competing orders at a fixed 
temperature T/t\\ = 0.4 in Fig. g](b). The FFLO state 
has a large value of entropy density due to the extra un- 
paired majority fermions, which interpolate between the 
BCS and the fully polarized normal state. This greatly 
increases the accessibility of the FFLO state in the cold 
atom optical lattices. The transition between the BCS 
and the FFLO states is first order as indicated by the 
discontinuity of entropy in Fig. |4] (b). 

At last, we discuss experiment realizations and detec- 
tions. The p-band fermion systems can be realized by 
first preparing enough number of atoms to fully fill the s- 
orbital band, thus the extra particles will fill the p-bands. 
The attractive interaction can be achieved through Fesh- 
bach resonances in lattices @, Eo| , whose strength can be 
tuned comparable to the band width of 4ty « 0.5^ at 
Vo/Er & 15 [30(| , but still small compared to band gaps 
which is around several Er. Our work predicts a large 
stable parameter regime for the FFLO states. These 
states can be detected by many methods [ill I231. 124} . such 
as the direct imaging of the density profile oscillations 
of each of the fermion components, the rf spectroscopy 
measurement on the collective modes, converting Cooper 
pairs into molecules and measuring their momenta, the 
shot-noise correlation of the Fermi momenta between k 
and — k ± Q, etc. In particular, the recent development 
of the in situ imaging methods with the single site res- 
olution [4]], EH can be used to accurately determine the 
spatial oscillation of the FFLO states. 

In summary, we have studied the competing orders 
among the FFLO, the BCS, and the normal states in 
the spin imbalanced p-orbital band systems in both 2D 
and 3D. The FFLO states are stabilized by the combined 
effects of the quasi- ID Fermi surfaces and the high di- 
mensional phase coherence built up by the inter-orbital 
interactions. The pairing density wave vectors are along 
the diagonal directions to facilitate the maximal inter- 
orbital pairing phase coherence. The FFLO states are 
robust with many realistic experimental effects including 
the confining trap, the small transverse 7r-bonding, and 
finite temperatures. It would be nice to realize the 2D 
and 3D stable FFLO phases in the p-orbital bands in op- 
tical lattices which have not been identified in solid state 
systems yet. 
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